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Abstract 

Novel detector schemes are proposed for the short-baseline neutrino experiments 
of next generation, aimed at exploring the large-Am^ domain of Vt oscilla- 

tions in the appearance mode. These schemes emphasize good spectrometry for 
charged particles and for electromagnetic showers and efficient reconstruction of 

_^ r^ry decays. The basic elements are a sequence of relatively thin emulsion 
targets, immersed in magnetic field and interspersed with electronic trackers, and a 
fine-grained electromagnetic calorimeter built of lead glass. These elements act as 
an integral whole in reconstructing the electromagnetic showers. This conceptual 
scheme shows good performance in identifying the r (quasi-)two-body decays by 
their characteristic kinematics and in selecting the electronic decays of the r. 
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1 Introduction 

Whether or not the leptonic number is strictly conserved is still an open problem of particle 
physics. The leptonic number may be violated, if the current neutrinos = e, fi, r) 

are not exactly massless, but rather represent linear combinations of the mass eigenstates 
Ui with masses rrii {i = 1, 2, 3): Ua = Uaii^i, where Uai is a unitary mixing matrix 
that is analogous to the Cabibbo-Kobayashi-Maskawa matrix for the quark sector. Then, 
as different mass components no longer propagate coherently, one current neutrino may 
oscillate in flight to another current neutrino ^ i^jj, a ^ p. Over a distance L 

from the point of emission, the transition probability is 

P{a ^ (3) = 5af3 - 4:'^UaiUpiUajUpjSin^{7rL/Xij), 

i>j 

where the oscillation length, Xij, is proportional to the neutrino energy E^y and is inversely 
proportional to the mass difference Amjj = mf — m'j : 

Note that in the considered case of just three Dirac neutrinos, we have only two indepen- 
dent mass differences Am^^-. Given the values of By and Am^, the transition probability 
reaches its first maximum at a distance L = 1.24 km (GeV) / Am^ (eV^) from the 
point of emission. 

Searches for neutrino oscillations are actively pursued in experiments with solar, at- 
mospheric, reactor, and accelerator neutrinos. These are characterized by very different 
L/Ey ratios and, therefore, are sensitive to very different values of Am^. A large deficit 
of muon neutrinos, as observed in atmospheric showers 0, suggests that they transform 
to another flavor with Am^ on the order of 10~^-10~^ eV^ and with an effective mixing, 
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defined as sin^(26') = 4:UaiU^iUajUf3j, between 0.5 and 1.0. On the other hand, an acceler- 
ator experiment using neutrinos from 7r+ and /i^ decays at rest has presented evidence 
for the transition with a much larger mass difference of Am^ ~ 1 eV^ and with 

a much smaller effective mixing of sin^(20) ~ 6 x 10^^. 

The former pattern of flavor-changing transitions between neutrinos, that is character- 
ized by a small Am^ and a large mixing, will be thoroughly investigated in long-baseline 
accelerator experiments like MINOS at Fermilab in which the muon neutrinos will 
travel a distance of 731 km before hitting the detector. The latter pattern of oscilla- 
tions, that involves Am^ ~ 1 eV^ and a small mixing, must be further studied in either 
the z/g and channels in short-baseline accelerator experiments that may 

accumulate sufficient statistics of neutrino interactions. 

Despite some appealing attempts to describe all available data as a whole we 
believe that the experimental situation is still too volatile to allow a coherent and reliable 
scheme of transitions among the three neutrino flavors. Still, there are sound qualitative 
reasons to believe that the transition driven by the bigger Am^^ must primarily manifest 
itself in the Vr channel. Qualitatively, one may expect that the mass hierarchy 

of neutrinos follows that of corresponding charged leptons, so that z/^ should be the 
most massive. Indeed, the GUT see-saw mechanism predicts that mv^lmy^jmy^ ~ 
ml^/rnl/m^. On the other hand, in analogy with the known pattern of quark mixings, 
one may naively expect that mixings between the neighboring neutrino flavors {i.e., z/g-z/^ 
and z/^-z/,-) should be the strongest, while that for z/g-z^r may be substantially weaker. 

Therefore, if the bigger Am^^- is on an order of 1 eV^ as suggested by the data of 
LSND 1^, one should look for the corresponding transition in the short-baseline 

(L ~ 1 km) and medium-baseline (L ~ 10 — 20 km) experiments at accelerators. Quite 
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importantly, the accelerators like SPS at CERN and Main Injector at Fermilab offer 
intense and almost pure beams of muon neutrinos with mean energies well above the 
threshold for r production, as required by experiments with Vj. appearance. 

The flux of r neutrinos generated via the flavor-changing transition — > Vj- can 
only be inferred from the number of r leptons produced in the charged-current reaction 
Vt-N T^X. The created r typically travels a distance of one millimeter in space 
(lifetime 3 x 10^^'^ s) and then decays either leptonically with emission of two neutrinos 
(r^ fi^uu, e^vv) or hadronically with emission of a single neutrino (r^ n^u, p^u, 
a5~z/, etc.). An efficient search for the z/^ Ur transition requires that the events of r 
production and decay be selected and analyzed on an event-by-event basis. Therefore, the 
production and decay vertices must be reliably resolved in the detector (whereby the short 
track of the r is reconstructed in space) and the decay particles must be identified. Of 
the available detector choices, only nuclear emulsion fully meets the above requirements 
for a massive target. 

The emulsion technique was pioneered by the Fermilab experiment E531 which, in 
its search for the z/^ transition p in the appearance mode, found no candidate 

events of r production and decay. For Am^ > 50 eV^, the mixing parameter sin^(2^) was 
restricted to be less than 5 x 10^1 At CERN SPS, the emulsion experiment CHORUS 



completed operation in 1997, but data processing is still in progress |Ty]. Provided that 
no candidate events are finally seen, they will tighten the E531 upper limit on sin^(2^) 
by another order of magnitude. 
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2 Short-Baseline Experiments of Next Generation 



Two short-baseline emulsion experiments of the next generation have been designed to 
achieve an order-of-magnitude increase in sensitivity over CHORUS and the electronic 
experiment NOMAD fill. At Fermilab, E803/COSMOS p|, |l|] was optimized for the 



intense neutrino beam of Main Injector - a 120-GeV proton accelerator to be commis- 



sioned in 1999. At CERN, TOSCA [|T|] will use the z/^, beam generated by the 350-GeV 
proton beam of the SPS accelerator. Like E531 and CHORUS, both are designed as 
hybrid detectors combining an emulsion target with a downstream electronic spectrome- 
ter. The candidates for r decays are directly observed in the emulsion target, while the 
spectrometer is used to select the events to be scanned in emulsion and to measure the 
momenta of secondary particles for kinematic analysis of r candidates. 

The baseline schemes of COSMOS and TOSCA are detailed in the respective proposals 



[0] and [|14[, and general characteristics of the two experiments are compared in Table 
1^. Either experiment will collect ten times more CC interactions of muon neutrinos than 
CHORUS, while rejecting the backgrounds to the r signal much more efficiently. The 
bigger target mass and higher beam energy of TOSCA are seen to be largely compensated 
by the higher repetition rate of Main Injector that will deliver to target some eight times 
more protons per year than CERN SPS. As a result, the two experiments may boast 
comparable sensitivities to the z/^ transition in the null limit, that is, assuming that 
no signals are finally detected. 

That the experiment be capable of convincingly interpreting and demonstrating even 
a relatively small r signal, should one show up, is perhaps much more important than just 
restricting the parameter space for the null hypothesis. For this, an appreciable fraction 



of the signal must be unambiguously reconstructed as r decays rather than the decays 
of anticharm produced by the antineutrino component of the beam. Such distinctive 
signatures can only be provided by the (quasi-) two-body semileptonic decays of the r : 

T- n-y {BR= 11.3%), 

T- ^ p-u, p- 7r"7r° (BR = 25.2%), 



Ti-n+n- {BR = 9.4%). 



In a (quasi-)two-body decay r~ /i~z/, "transverse mass" is defined as Mt = \J'm\ + Pr+ 
PTi where mt and px are the h~ mass and transverse momentum with respect to r 
direction. The two-body kinematics dictate that the Mt distribution should reveal a 
very distinctive peak just below M^"^ = m,-. This "Jacobian cusp", that may provide 
a characteristic signature of the r, rapidly degrades with increasing /S.p/p for the decay 
products (see [|12|, |1^ and the figures below). Needless to say, analyzing the r~ p~u 



decays also requires good reconstruction of 7r° 77 decays in the detector. The Mt 
technique for identifying massive parents by two-body decays in emulsion was proven by 
observing the relatively rare decays D^(1968) — > p^v against a heavy background from 
other decays of charm |15|. 



In its baseline form that emphasizes good spectrometry for charged particles {Ap/p ~ 
0.03 disregarding the error from multiple scattering in emulsion), COSMOS will be able 
to observe a distinct Jacobian cusp in the Mt distribution for the pionic decay 
TT~v. The decay r~ p~v will be analyzed less efficiently because of difficulties in 
reconstructing the vr'' 77 decays in thick emulsion (three radiation lengths in the 
baseline design). Thick emulsion will also complicate the analysis of electronic decays 

T- e'uu {BR = 17.8%) 



and subsequent comparison with the muonic decays 

r- uu (Si? = 17.4%). 

In TOSCA, the analysis of (quasi-)two-particle decays of the r is severely hampered by 
inferior spectrometry for charged particles (Ap/p ~ 0.1) and by the fact that the decays 
7r° 77 are not reconstructed at all. The one-prong semileptonic decays 

T- 7r-(n7r°)i/ (BR ~ 50%) 

are treated only inclusively, i.e., without selecting the individual (quasi-)two-particle chan- 
nels. In this sense, in its baseline form TOSCA is a "counting experiment" that only com- 
pares the number of observed kinks with predicted backgrounds from anticharm decays 
and from pion interactions in emulsion. Thus, despite a somewhat higher sensitivity to 
the i/^ — > i^T transition in the null limit, TOSCA may have less discovery potential than 
COSMOS. 

3 Alternative Schemes for COSMOS and TOSCA 

In this paper, we propose alternative conceptual schemes for COSMOS and TOSCA, 
aimed at enhancing the discovery potentials of both experiments through improving the 
analysis of (quasi-)two-body decays of the r. Accordingly, the proposed schemes empha- 
size good spectrometry for charged particles and electromagnetic showers and efficient 
reconstruction of 7r° — > 77 decays. 

Either scheme involves a sequence of relatively thin emulsion targets, immersed in 
magnetic field and interspersed with electronic trackers, and a fine-grained electromag- 
netic calorimeter (EMCal) downstream of the targets. As the total thickness of emulsion 
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amounts to several radiation lengths, most photons originating from the targets will pro- 
duce electromagnetic showers. The adopted strategy is to reassemble the parent photon 
from unconverted daughter gammas (detected in the EMCal) and from conversion elec- 
trons that are momentum-analyzed in the instrumented gaps. By additionally sampling 
the shower in between the targets that act as converters, we are able to collect a signifi- 
cant fraction of parent energy despite relatively large total thickness of emulsion. Thus, 
in reconstructing electromagnetic showers the targets, trackers, and the EMCal act as 
an integral whole. This strategy develops the original philosophy of COSMOS and is 
very alien to the baseline design of TOSCA in which the six target modules are virtually 



independent of each other 1 14 ] 



The strategy of assembling the parent photon from "pieces" dictates that the EMCal 
should show good performance at energies down to some 100 MeV in a high-multiplicity 
environment (10-20 showers per event). Of the existing options, these requirements are 
best met by the technique of Cherenkov lead glass. For either COSMOS and TOSCA, 
we assume a round-shaped EMCal with 3.5-m diameter, built of lead-glass cells with 
transverse granularity of 4.25 by 4.25 cm and with thickness of 15 radiation lengths. 
The cells are read out by photomultipliers. We have ample experience in constructing 
and operating very similar devices, and the response of the proposed EMCal is fully 
understood through extensive simulations and calibrations [0. It is characterized by an 



energy resolution of 6E (GeV) = 0.05v-E' + 0.02ii^, while the position resolution in either 
coordinate is close to 2 mm. 

In the proposed schemes, as in the baseline design of COSMOS, electronic track- 
ing largely relies on the multisampling (jet) drift chambers that provide excellent pat- 
tern recognition, redundancy, and two-track resolution. Following the original design of 
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TOSCA, we also include two planes of silicon-microstrip detectors, providing two inde- 
pendent XY measurements to define a track segment, just downstream of each target (at 
10 and 50 mm from bulk emulsion, respectively). These are primarily aimed at facilitat- 
ing the extrapolation of tracks into bulk emulsion [|T7|, but also affect the spectrometry 
of these tracks. In the simulations described in subsequent sections, we assume that the 
drift chambers and silicon-microstrip planes have spatial resolutions of 150 and 20 fim, 
respectively. 

In either scheme detailed below, we assume that the targets and trackers are immersed 
in uniform magnetic field of 0.5 Tesla that is perpendicular to beam direction. The 
geometry of the magnet is not specified. (Note however that TOSCA will use the existing 
dipole magnet with useful volume of 3.5 x 3.5 x 7.5 m^ and field of up to 0.7 Tesla, that 
fully meets our requirements. At the same time, our scheme for COSMOS will require a 
bigger magnet than that foreseen in the original proposal.) Likewise, we do not specify 
the geometry of the muon identifier that is not relevant to our analysis. 

The proposed configuration of COSMOS is shown in Fig. [l|. The three emulsion 
targets with areas of 180 x 180 cm^ and thicknesses of 3 cm (~ one radiation length) 
are located at 2; = 0, 120, and 240 cm along the beam direction. The front face of the 
EMCal is at 2 = 540 cm, or far enough from the targets for efficiently reconstructing the 
71^ — > decays. Also shown are the pairs of silicon-microstrip planes just downstream of 
each target and the drift chambers with active area of 240 x 240 cm^ that instrument the 
120-cm-wide gaps between the targets and the 300-cm-wide gap between the downstream 
target and the EMCal. The total thickness of emulsion is the same as in the baseline 
scheme (9 cm), but now we are able to sample electromagnetic showers every radiation 
length of emulsion. The total amount of emulsion is 1112 kg, or some 30% more than in 
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the baseline design. 

For TOSCA, we propose a very similar configuration that is also depicted in Fig. |1]. 
Here, the number of targets is increased to four (at z = 0, 150, 300, and 450 cm) and 
their thicknesses are increased to 4 cm, thus bringing the total thickness of emulsion to 
over five radiation lengths. Compared to COSMOS that will operate at a lower beam 
energy (see Table |l]), the targets are driven farther apart to provide a longer lever arm for 
momentum analysis. The targets and the drift chambers have the same areas as in the 
COSMOS configuration (180 x 180 and 240 x 240 cm^, respectively). That the EMCal is 
now placed farther downstream of the last target (front face at z = 950 cm) is dictated 
by higher mean energy of the 7r° mesons to be reconstructed by TOSCA. Compared to 
the baseline design of TOSCA, the total amount of emulsion is less by some 17%, while 
the total area of silicon-microstrip planes is virtually the same. 



4 Simulating r Signatures in the Detector 

For either scheme detailed in the previous section, we have performed GEANT-based 
simulations of detector response to the process of r production and decay through various 
channels. In our calculations, we assume that the energy spectra of incident r neutrinos 
are exactly proportional to those of original z/^ beams from the respective accelerators 
{i.e., from Main Injector for COSMOS and from CERN SPS tuned to 350-GeV proton 
energy for TOSCA, see Table p. This corresponds to u^^ — > oscillations driven by a 
large mass difference (Am^ > 30 and 50 eV^ for COSMOS and TOSCA, respectively). 
The radial distribution of either beam at detector position is taken into account. The 
computed energy spectra of i/^-induced CC interactions for either experiment, that are 
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affected by the energy threshold for r production, are shown in Fig. ^ 

Note that a reahstic simulation should include the experimental restrictions dictated 
by the emulsion technique: 



• Because of Coulomb scattering, soft and broad-angle tracks can neither be recon- 
structed in emulsion nor be linked to the tracker. Therefore, a charged daughter of 
the T can only be detected and momentum-analyzed, if its lab. momentum exceeds 
some 1 GeV and if its emission angle is within some 400 mrad. 

• In a one-prong decay (like l^uu, 7r~z/, or P~^)i the angle between 
the r and charged-daughter directions must exceed some 10 mrad for the "kink" on 
the r track to be detected in emulsion. 

• At low transverse momenta, the bulk of detected kinks are due to decays of strange 
particles and to pion scatterings on the nuclei of emulsion without visible breakup 
of the nucleus (the so-called "white stars"). These background kinks are largely 
removed by the cut pt > 250 MeV on transverse momentum of the charged daughter 
with respect to parent-r direction. 

A major source of background to the r signal is anticharm production in those 
antineutrino-induced CC events in which the primary charged lepton (either a or 
a e"*") has been misidentified in the detector. That the r, unlike anticharm, is emitted 
back-to-back with primary hadrons in the transverse plane allows to suppress the an- 
ticharm background by additional kinematic cuts (see |l^, |14|) that are not discussed 
in this paper. 
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4.1 Simulation Procedure 



The first stage is to generate the CC coUisions of r neutrinos with free nucleons, Vt-N 
T~X, with subsequent leptonic and semileptonic (quasi)-two-particle decays of the emitted 
r: T~ — > l~i'P, 7r~i', p~v, and a^v. Both the quasielastic and deep-inelastic processes of 
r production are simulated. For quasielastics, the hadronic tensor is expressed through 
the known formfactors for the weak n — > p transition. (Excitation of baryonic resonances 
like A"*" and A"*""*" is not simulated explicitly, but is accounted for by simply scaling up by 
a factor of 2.2 the cross section for Vt-N — > t'N' on an isospin-averaged nucleon.) For 
inelastic production of the r, we rely on the naive parton model. 

The process of r production and decay is treated as a whole using the narrow-width 
approximation: for each decay channel, the cross section is computed from an overall 
Feynman graph with two weak vertices in which 5{q^ — m^) is substituted for the r 
propagator. In this way, we are able to account for polarization of the r that strongly 
affects the angular distributions of secondary particles in the r frame: thus, in the decay 
T~ — > 7r~i/ the 7r~" is not isotropic in r frame but rather travels in the backward direction. 
(Note that polarization of the secondary resonances, p~ and aj^, is not taken into account, 
i.e., tertiary pions are assumed to be isotropic in the rest frame of the parent resonance.) 
The lineshapes of intermediate resonances are generated in Gaussian forms. 

For a given decay channel of the r, the output of the generator is a sequence of 
production-decay topologies (including tertiary pions from decays of secondary reso- 
nances) accompanied by weights that represent the computed matrix elements. Cor- 
respondingly, all distributions are plotted using these weights. At the moment, the indi- 
vidual hadrons of the primary hadron jet are not generated, and are treated as undetected. 
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Thus generated r events are then propagated through the detector using GEANT. 
Rather than fully digitize the response of electronic trackers (drift chambers and silicon 
microstrips) , the individual hits are randomly distributed around central positions ac- 
cording to corresponding resolutions. Then, particle momentum is obtained by fitting 
the "smeared" trajectory in magnetic field. Likewise, the hits in individual cells of the 
EMCal are not fully simulated at this early stage. Instead, the energy and position of a 
shower are smeared according to aforementioned resolutions. The problem of overlapping 
showers is not yet addressed. 

For a given decay channel of the r, all histograms are normalized to the total number 
of occurring decays (so that the net contents are just the accepted fraction of all r decays 
in this channel). With this convention, the acceptance of a selection becomes immediately 
obvious without further normalization. 

4.2 Detecting the Decay r~ ^ tt^u 

The momentum resolution for the tt^ with > 1 GeV from the decay —>■ n^u, as 
illustrated in Fig. ^, on average amounts to some 2.7% for either detector. The unsmeared 
transverse mass, Mt = + pf^ + px, is plotted for all generated r~ 7r~u events in 

Fig. |. Also shown are the corresponding smeared distributions for those decay pions 
that pass the procedure-driven "default" selections in lab. momentum (p^ > 1 GeV), in 
emission angle relative to incident beam (0 < 400 mrad), in kink angle {6rn > 10 mrad), 
and in transverse momentum to r direction {pt > 250 MeV). The Jacobian cusp of the 
original Mt distribution is but slightly degraded by apparatus smearings of the proposed 
COSMOS and TOSCA detectors. Therefore, in both detectors the produced r leptons 
will be efficiently tagged by tt^z/ decays. 
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4.3 Detecting the Decay r ^ z/, tt" ^ 77 

In this subsection that deals with the decay —>■ p^u, the "default" selections for the 
7r~ from ti^tt^ (namely, p^r- > 1 GeV, > 400 mrad, 6^^^ > 10 mrad, and p'^ > 250 

MeV) are implicitly included in all distributions and quoted acceptances. For the tt", we 
use the values of emission angles as measured near the production point in emulsion. 

We assume that the conversion point of the photon from tt^ 77 is found and 
measured in emulsion and, therefore, its direction is precisely known. The energy of 
the parent photon is reconstructed by adding up the energies of unconverted daughter 
gammas (by hits in the EMCal) and of conversion electrons that are momentum-analyzed 
in the tracker. 

In this scheme of reconstruction, there is a potential danger of double counting: a 
conversion electron, that has already been momentum-analyzed in the gap between the 
two targets, may shower in the next target producing further conversion electrons and 
gammas that should no longer be counted. We assume that, given a relatively large gap 
between successive targets, these "unwanted daughters" can be distinguished by their 
large displacement from the shower axis. Whether or not this can be done experimentally 
is currently under investigation, and our preliminary results are very encouraging. But 
for the moment, we realize this on GEANT level by "stopping" conversion electrons just 
before they hit the next target. 

By either allowing or excluding double counting, we obtain two different estimates of 
the parent-photon energy. For either approach, the ratio E'^'^'^/ between the estimated 
and true values of energy is plotted in Fig. ^ for photons originating from different 
emulsion targets of reconfigured COSMOS. Allowing "double counting" in is seen to 
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result in overestimating the true energy for a significant fraction of photons from vr" ^ 77. 
Again for either estimate of -E"*, plotted in Fig. ^ is the visible two-photon mass m^^ for 
events occurring in different targets of COSMOS. Quite predictably, "double counting" in 
j^ms (^QgYQ^dQg the 7r° signals in the two upstream targets. The very similar distributions 
of q£ different targets of reconfigured TOSCA are not illustrated. 

With "double counting" excluded, the overall m^^ distributions for the reconfigured 
COSMOS and TOSCA detectors are shown in Fig. |^. Fixing the mass window for 7r° 
selection requires detailed simulations of detector response to primary hadrons and to 
background processes, and therefore is impossible at this stage. For several realistic mass 
windows for m^^ (135 ± 25, 20, and 15 MeV), the acceptances of reconfigured detectors 
to —>■ p^v decays are listed in Table |^ together with those of the original COSMOS 
detector |13[. (Note that reconstruction of 7r° 77 decays, that is essential for selecting 
the T~ p'u decays, is virtually impossible in the baseline design of TOSCA |T^.) 
Either of the proposed detectors is seen to have a substantially higher acceptance to 
T~ p~v than the baseline detector of COSMOS. 

Tentatively selecting the 7r° 77 decays in the m^^ window of 135 ± 25 MeV, we 
then plot the mass of the 7r~7r° system, as reconstructed in either detector (see Fig. 
1^). And finally, in Fig. | we plot the reconstructed transverse mass for — > p^z/, 
Mt = \J + + Pt- As the Jacobian structure of the original Mt distribution is not 
destroyed by apparatus smearings, we may conclude that both proposed detectors will 
efficiently tag r leptons by reconstructing the ^ p^u decays. 
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4.4 Detecting the Decay r ai ly, ai ^ tt tt+tt 

This decay proceeds through an intermediate p° (a|f —>■ p^n^) and has a branching fraction 
of 9.4%. The three-prong decay channel is plagued by the background from three-prong 
collisions of secondary hadrons with the nuclei of emulsion, including coherent production 
of the ai state by pions. Still, it may be worthwhile to check whether or not the expected 
Mt distribution is distinctive. 

Using the procedure-driven selections of > 1 GeV and < 400 mrad for all three 
pions and determining pion directions near the production point in emulsion, in Fig. ^ 
we plot the reconstructed transverse mass for the decay o-i^, o^i tt^tt+tt^ : 

Mt = \J + p\ + pt- (Note that no Or-K selection is required for observing a three- 
prong decay in emulsion.) That the very distinct Jacobian structure of the original Mt 
distribution persists in the smeared plots is very encouraging: backgrounds from pion 
dissociation and from coherent ai production should be amassed at much lower values of 
Mt- Whether or not the proposed detectors are good enough for separating the decays 

^ can only be decided by detailed comparisons with the simulated three-prong 
background. 

4.5 Detecting the decay r e^vv 

Of the decays that are not in the (quasi-) two-body category, we consider the decay 
e^vv that, compared to pTvv^ poses experimental problems due to electrons 

showering in thick emulsion. Owing to thinner emulsion targets and to good reconstruc- 
tion of shower energy, we may detect the r~ e^vv decays more efficiently than the 
baseline apparata of COSMOS and TOSCA. 
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In reconstructing the electron from r~ e~viJ, we follow a strategy analogous to 
the one used for photons: the energy is reassembled from pieces, using magnetic analysis 
for all positive and negative electrons leaving the emulsion and hits in the EMCal for 
the unconverted daughter gammas. Thus reconstructed energy is denoted as E'l'^'^ . Then, 
we are able to estimate the transverse momentum with respect to r direction as p^** = 
i?g** sin6'T-e, where Ore is the angle of the observed kink. 

That the sign of the generic electron be correctly assigned is important for rejecting 
the background from semileptonic decays of charm produced by neutrinos. Therefore, the 
e~ track must be followed in emulsion down to the exit and reliably linked to electronic 
trackers. For this to be possible, we assume that the energy of the generic-e~ track at 
the exit point from the target of origin, Ef^^lj., must exceed 1 GeV. In a thick target, the 
generic e~ may fail this requirement by losing too much energy through bremsstrahlung, 
even though the total energy of the shower can be estimated. 

Thus, a candidate event for — > e~z/i^ must satisfy the following conditions: Ef^^^f. > 
1 GeV, 6e < 400 mrad, pp"" > 250 MeV, and Ore > 10 mrad. For the r" e~vv events 
that survive these minimum selections in the detector, we plot the ratio between the 
reconstructed and original energies of the e~, in pig. ^. The reconstruction 

of electron energy is good in both detectors proposed. For either the reconfigured and 
baseline designs of COSMOS and TOSCA, the fractions of all generated — ^ e^vu 
events that survive the above selections are listed in Table ^. The reconfigured detectors 
are seen to have higher acceptances to the decay e^vv than the baseline schemes 

for COSMOS and TOSCA. 
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5 Summary 



The large-Am^ domain of — > v^. oscillations may be further probed by the next gener- 
ation of short-baseline neutrino experiments using the emulsion technique for detecting r 
decays: COSMOS at Fermilab and TOSCA at CERN. In this paper, we propose alterna- 
tive schemes of these experiments that emphasize good spectrometry for charged particles 
and for electromagnetic showers and efficient reconstruction of 7r° — > 77 decays. Either 
configuration features a sequence of relatively thin emulsion targets, immersed in mag- 
netic field and interspersed with electronic trackers, and a fine-grained electromagnetic 
calorimeter built of lead glass. These elements act as an integral whole in reconstructing 
the electromagnetic showers. This conceptual scheme shows superior performance in iden- 
tifying the r (quasi-)two-body decays by their characteristic kinematics and in selecting 
the electronic decays of the r. 



19 



References 

[1] B. Pontecorvo, Zh. Eksp. Teor. Fiz. 7, 172 (1958); Sov. Phys. JETP 6, 429 (1958). 

[2] M. C. Gonzales-Garcia, H. Nunokawa, O. L. G. Peres, T. Stanev, and J. W. F. Valle, 
Update on Atmospheric Neutrinos, Phys. Rev. D58, 033004 (1998). 

[3] C. Athanassopoulos et al, LSND Coll., Phys. Rev. Lett. 77, 3082 (1996). 

[4] E. Abies et al, MINOS Coll., Fermilab Proposal P875 (1995). 

[5] P.F. Harrison, D.H. Perkins, and W.G. Scott, Phys. Lett. B349, 137 (1995). 

[6] A. Acker and S. Pakvasa, preprint hcp-ph//9611423 (1996). 

[7] M. Gell-Mann, P. Ramond, and R. Slansky, in Supergravity, eds. P. van Nieuwen- 
huizen and D. Freedman (North Holland, Amsterdam, 1979), p. 315. 

[8] N. Ushida et al, E531 Coll., Phys. Rev. Lett. 57, 2897 (1986). 

[9] M. de Jong et al, CHORUS Coll., preprint CERN-PPE/93-131 (1993); E. Eskut et 
al, CHORUS Coll., Nucl. Instr. and Meth. A401, 7 (1997). 

[10] E. Eskut et al, et al, CHORUS Coll., Phys. Lett. B424, 202 (1998); E. Eskut et al, 
CHORUS Coll., preprint CERN-PPE/98-73 (1998). 

[11] J. Altegoer et al, NOMAD Coll., Nucl. Instr. and Meth. A404, 96 (1988); J. Altegoer 
et al, NOMAD Coll., preprint CERN-EP/98-57 (1998). 

[12] K. Kodama et al, E803/COSMOS Coll., Fermilab Proposal P803 (1993). 

[13] K. Kodama et al, E803/COSMOS Coll., Update Report on Fermilab E803/COSMOS 
(1995). 

20 



[14] A.S. Ayan et al, TOSCA Coll., Letter of Intent, CERN-SPSC/97-5, SPSC/I213 
(March 1997). 

[15] K. Kodama et ai, E653 Coll., preprint DPNU-96-33 (June 1996). 

[16] R. Rameika, A. Asratyan, M. Balatz, G. Davidenko, A. Dolgolenko, G. Dzyubenko, 
V. Kaftanov, M. Kubantsev, A. Sitnikov, V. Verebryusov, and V. Vishnyakov, Lead- 
Glass Option for the COSMOS Electromagnetic Calorimeter, E803/COSMOS inter- 
nal report (February 1996). 

[17] S. Aoki et al., Nucl. Instr. and Meth. B51, 466 (1990). 



21 



Quantity Compared 


COSMOS 


TOSCA 


Proton energy, GeV 


120 


350 


Protons per cycle 


4 ■ 10^^ 


6 ■ 10^^ 


Cycle time, s 


1.9 


19.2 


Protons per year 


3.7-1020 


4.3 ■ 10^9 


Length of decay channel, m 


800 


414 


Distance from target to detector, m 


960 


806 


Emulsion mass, kg 


865 


2784 


Tn 1 • 9 

iimuision area, m 


1.80 X 1.40 


1.44 X 1.44 


Emulsion thickness, cm 


2 X 4.5 


6 X 6.0 


{E) of z/^ beam (by flux), GeV 


12 


25 


{E) of z/^-induced CC events, GeV 


18 


37 


{E) of z/T--induced CC events, GeV 


24 


54 


{L/E), km/GeV 


0.023 


0.012 




0.27 


0.50 


S (BR X efficiency) 


0.114 


0.104 


Total number of z/^-induced CC events 


8.1 ■ 10^ 


6.0 ■ 10^ 


Total number of background events 


1.3 


1.3 


Total background events, 90% CL 


3.5 


3.5 


Reach in sm'^{29) for large Am^ 


2.8 ■ 10"^ 


1.5 ■ 10"^ 


Reach in Am^ for maximum mixing, eV^ 


0.10 


0.10 



Table 1: General parameters of the baseline designs of COSMOS ||T3l and TOSCA [jl4| 
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Mass window for 


Our design 


Our design 


Baseline design 


selecting the 7r° 


for COSMOS 


for TOSCA 


for COSMOS 


Unrestricted 


0.40 


0.42 


0.31 


135 ± 25 MeV 


0.31 


0.28 


0.19 


135 ± 20 MeV 


0.29 


0.25 


0.18 


135 ± 15 MeV 


0.26 


0.21 


0.16 



Table 2: The fractions of all r — > p i/, p — > tt 7r° decays that survive different tt' 
selections (as indicated) and the "default" selections for the 7r~. 
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Our design 
for COSMOS 


Our design 
for TOSCA 


Baseline design 
for COSMOS 


Baseline design 
for TOSCA 


0.55 


0.56 


0.31 


0.42 



Table 3: The fractions of all r — > e vv decays that survive the minimum selections in 
the detector, as detailed in the text. 
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Figure Captions 

Fig. |T]. The configurations proposed for COSMOS (top) and TOSCA (bottom). 
Overlaid are tlie p^v decays in eittier detector. 

Fig. ^. Tlie computed energy spectra of ^'T--induced CC events in COSMOS (a) and 
TOSCA (b). The absolute normalization is arbitrary. 

Fig. ^. In the decay r~ ti~v, the ratio between the measured and true values of vr" 
momentum, R = for the proposed schemes of COSMOS (a) and TOSCA (b). 

Here and in all subsequent figures, the smearing due to Coulomb scattering in emulsion 
is fully taken into account. 

Fig. ^. Transverse mass Mt = \Jpt + 'mfi + Pt for the (quasi-)two-body decays 
—>■ h^v with h~ = tt^ (left-hand column), = ^ n^n^ (middle column), and 
h~ = ai ^ 7i~7T~^'K~ (right-hand column). The unsmeared My distributions for all 
generated events in each channel are shown in the top row. The smeared distributions for 
events detected by the reconfigured COSMOS and TOSCA are shown in the middle and 
bottom rows, respectively. 

Fig. ^. The ratio R = EJ!^^'^ / Et^'"^'^ between the estimated and true values of energy 
for photons originating from different emulsion targets of reconfigured COSMOS, with 
"double counting" in E'™* either allowed (the top row) or forbidden (the bottom row). 
The three columns stand for the three targets of the proposed detector. 

Fig. ||. The reconstructed two-photon mass m^^ for the 7r° — > 77 decays occurring 
in different emulsion targets of reconfigured COSMOS, with "double counting" in 
either allowed (the top row) or forbidden (the bottom row). The three columns stand for 
the three targets of the proposed detector. 
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Fig. 1^. The masses m^^ (top row) and m^rTr (bottom row), as reconstructed in the 
reconfigured COSMOS and TOSCA detectors (the left-hand and right-hand columns, 
respectively) excluding double counting in i?^**. Note that the lineshape of the original 
p~ has been generated as a Gaussian with cr = 75 MeV. 

Fig. 1^. The ratio between the reconstructed and original energies of the e~ from 
r- ^ e-i/z/, R = El^'/El^''^, for the proposed configurations of COSMOS (a) and TOSCA 
(b). 
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Figure 2: The computed energy spectra of z/T--induced CC events in COSMOS (a) and 
TOSCA (b). The absolute normahzation is arbitrary. 
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Figure 3: In the decay t~ tt^v, the ratio between the measured and true values of tt" 
momentum, R — p'^^^ / p*^'^^ , for the proposed schemes of COSMOS (a) and TOSCA (b). 

Here and in all subsequent figures, the smearing due to Coulomb scattering in emulsion 



is fully taken into account. 
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Figure 4: Transverse mass Mt = sjpr + + Pt for the (quasi-)two-body decays t~ — >• 
h'u with h~ — Tr~ (left-hand column), h~ — p~ ^ 7r~7r° (middle column), and h~ — 
ai — > 7r~7r"'"7r~ (right-hand column). The unsmeared Mt distributions for all generated 
events in each channel are shown in the top row. The smeared distributions for events 
detected by the reconfigured COSMOS and TOSCA are shown in the middle and bottom 
rows, respectively. on 



TARGET 1 



TARGET 2 



TARGET 3 




MEAN 0.96 
RMS 0.13 
ALLCHAN 0.30 




0.07 i- 

0.06 j- 

0.05 !■ 

0.04 j- 

0.03 j- 

0.02 j- 

0.01 j- 

I- 



MEAN 0.85 
RMS 0.19 
ALLCHAN 0.24 




0.07 i- 

0.06 j- 

0.05 j- 

0.04 j- 

0.03 j- 

0.02 j- 

0.01 j- 

.- 



MEAN 0.90 
RMS 0.17 
ALLCHAN 0.28 




0.07 i- 

0.06 j- 

0.05 j- 

0.04 j- 

0.03 j- 

0.02 j- 

0.01 j- 




MEAN 0.97 
RMS 0.12 
ALLCHAN 0.30 




Figure 5: The ratio R — E^'^^ / E^'^^ between the estimated and true values of energy 
for photons originating from different emulsion targets of reconfigured COSMOS, with 
"double counting" in i^""* either allowed (the top row) or forbidden (the bottom row). 
The three columns stand for the three targets of the proposed detector. 
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Figure 6: The reconstructed two-photon mass rrij^ for the tt 77 decays occurring in 
different emulsion targets of reconfigured COSMOS, with "double counting" in either 
allowed (the top row) or forbidden (the bottom row). The three columns stand for the 
three targets of the proposed detector. 
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Figure 7: The masses m-y-y (top row) and mTTTr (bottom row), as reconstructed in the 
reconfigured COSMOS and TOSCA detectors (the left-hand and right-hand columns, 
respectively) excluding double counting in E!^^^. Note that the lineshape of the original 
p" has been generated as a Gaussian with a — 75 MeV. 
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TT u, the ratio between the measured and true values of tt 
EkgrniKHieaiii^ ft^(^hfp^''^ , for the proposed schemes of COSMOS (a) and TOSCA (b). 

Here and in all subsequent figures, the smearing due to Coulomb scattering in emulsion 



is fully taken into account. 
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Figure 4: Transverse mass Mr = \JPt + + Pt for the (quasi-)two-body decays t~ — >• 
h~v with /i~ = 7r~ (left-hand column), h" — p~ ^ 7r~7r° (middle column), and = 
of — > TT^TT+TT" (right-hand column). The unsmeared Mt distributions for all generated 
events in each channel are shown in the top row. The smeared distributions for events 
detected by the reconfigured COSMOS and TOSCA are shown in the middle and bottom 
rows, respectively. qq 
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Figure 5: The ratio R — E^'^^ / E^'^^ between the estimated and true values of energy 
for photons originating from different emulsion targets of reconfigured COSMOS, with 
"double counting" in E"""* either allowed (the top row) or forbidden (the bottom row). 
The three columns stand for the three targets of the proposed detector. 
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Figure 6: The reconstructed two-photon mass rrij^ for the tt 77 decays occurring in 
different emulsion targets of reconfigured COSMOS, with "double counting" in either 
allowed (the top row) or forbidden (the bottom row). The three columns stand for the 
three targets of the proposed detector. 
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Figure 7: The masses m-y-y (top row) and mTTTr (bottom row), as reconstructed in the 
reconfigured COSMOS and TOSCA detectors (the left-hand and right-hand columns, 
respectively) excluding double counting in Note that the lineshape of the original 

p" has been generated as a Gaussian with a — 75 MeV. 
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Figure 8: The ratio between the reconstructed and original energies of the e from r — > 
e-uu, R = Ef'/E^;""^, for the proposed configurations of COSMOS (a) and TOSCA (b). 
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